From Chiral quark dynamics with Polyakov loop 
to the hadron resonance gas model 1 2 

E. Ruiz Arriola *, E. Megias t and L. L. Salcedo* 

* Departamento de Fisica Atomica, Molecular y Nuclear and Instituto Carlos I de Fisica Teorica y 

Computacional. Universidad de Granada, E-18071 Granada, Spain 
^Grup de Fisica Teorica and IFAE, Departament de Fisica, Universitat Autdnoma de Barcelona, 

Bellaterra E-08193 Barcelona, Spain 

Abstract. Chiral quark models with Polyakov loop at finite temperature have been often used to 
describe the phase transition. We show how the transition to a hadron resonance gas is realized 
based on the quantum and local nature of the Polyakov loop. 
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INTRODUCTION 

The cross-over between the chiral symmetry restoration and deconfinement at a common 
critical temperature T c ~ 200MeV was a first 01 0, 0] and by now firmly established 
prediction of lattice QCD [0, 0]. The transition point is characterized by two order 
parameters. On the one hand the quark condensate (qq) vanishes smoothly. On the 
other hand the Polyakov Lj = (tr c e' A °/ r ) /N c loop, where Aq is a gluon field in the 
adjoint representation of the SU(iV c ) gauge group, jumps smoothly from zero to one 
signaling the breaking of the center symmetry 7*(N C ) M\. Of course, (qq) and Lj are 
true order parameters in the opposite limits of vanishing quark masses (chiral limit) 
and for infinitely quark masses (gluodynamics) respectively, while the true cross-over is 
defined by inflexion points of both (qq) and Lj. The expectation that a phase transition 
between a hadronic phase to the quark-gluon plasma phase could be observed on the 
laboratory has inspired a wealth of work in recent years 170 . 

There have been many attempts to model the chiral-deconfinement cross-over, the 
main difficulty lies in properly combining the relevant degrees of freedom for the 
corresponding order parameter; while well below T c hadrons provide a complete basis 
of states, much above T c just quarks and gluons seem the adequate basis. The cross-over 
region seems difficult as it marks the coexistence of both degrees of freedom. In this 
contribution we will focus on the low temperature region of chiral quark models where 
the crossover is known to occur at higher temperature. 
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THE HADRON RESONANCE GAS MODEL 



On general grounds, and because quarks and gluons are confined, one expects that 
at sufficiently low temperatures any observable in QCD should be represented by the 
relevant hadronic colour neutral states. 

For the vacuum energy at finite temperature, and more specifically the trace of the 
energy momentum tensor, A(T) = (e — 3p)/T 4 , the results found on the lattice are 
expected to be represented by an interacting gas of low-lying stable hadrons (for light u 
and d quarks it would just be a gas of it, N and N states). Most of the interactions in the 
scattering region generate resonances (such as p, ft), A, etc.) which could be represented 
as narrow states provided the ratio T/M is small. This is consistent with the large N c 
limit expectation that T/M = @{N~ l ) while experimentally T/M = 0.12(8} M,m for 
both mesons and baryons listed in the Particle Data Group (PDG) booklet [10]. Thus, in 
the Hadron Resonance Gas (HRG) model the interactions are represented by a bunch of 
narrow-looking resonances whose partition function is given by flll fl2l. fl3lfl4l. l5ll. 



logZ=-/^£W„log(l-f„ e -^/-), (1) 

with g a = (2J a + l)(2r a + 1) the degeneracy factor, £ a = ±1 for bosons and fermions 
respectively, M a the hadron mass and z a the fugacity. The PDG states ifioll saturate 
lattice calculations [fl5lfl6ll Fl as also found within a strong coupling expansion for heavy 
quarks in lll7n . 

This HRG representation is less obvious for QCD operators involving just gluon 



fields. However, we have recently shown I118h that a hadronic representation of the 
Polyakov loop is given by 

where g a are the degeneracies and A a are the masses of hadrons with exactly one heavy 
quark (the mass of the heavy quark itself being subtracted). The comparison with the 
spectrum with u,d,s light quarks and one extra heavy quark turns out to be rather 
satisfactory and fairly independent on taking charm, bottom or truely infinite heavy 
quarks. It is also intriguing since these calculations might provide a handle on deciding 



the existence of exotic multiquark states [18] 



These HRG approximations are expected to hold at sufficiently low temperatures and 
agreement with lattice data is observed within the finite lattice uncertainties. In any 
case, it turns out that many states are needed to saturate both the trace anomaly and the 
Polyakov loop at temperatures below T c as there are no significant gaps in the spectrum. 



3 A temperature shift of about T = 10 - 20MeV is required; AhrgC^ - T Q ) = Aqcd(T) in Ref. |QJ]. 



POLYAKOV-CHIRAL-QUARK MODELS AND HEAVY QUARKS 



An effective and phenomenologically successful approach to the physics of the phase 
transition is provided by chiral quark models coupled to g luon fields in the form of 



27, 28 



3211 . Because 



a Polyakov loop |fl9l |2Cl |2ll IH |23l, [24 
most often the venerable Nambu-Jona-Lasinio model has been used, these models are 
referred to as PNJL models. In this contribution we discuss how PNJL models may be 
represented as a HRG at low T. Many works remain within a mean field approximation 
and assume a global Polyakov loop. As we have repeatedly criticized in our previous 
works [SII EM IH this raises the theoretical problem of the undesirable ambiguity 
of group coordinates on the one hand, but also the practical problem that the adjoint- 
representation provides a non-vanishing value for the Polyakov loop, contradicting 
lattice simulations. These difficulties may be overcome ||2ll. l22 |l by recognizing the local 
and quantum nature of the Polyakov loop. 



We start out from the partition function motivated in [211, 122L |33i [34J] 



Z = jDae- s ^ a l 



(3) 



where £1 = e iA °/ r and DQ is the invariant SU(iV c ) Haar group integration measure, for 
each SU(A^ C ) variable Q.(x) at each point .v. Here the action is 



S(T,Q.)=S q (T,0)+S G (T,a.). 



(4) 



The fermionic contribution depends on the quarks (and anti-quarks) is obtained from the 
corresponding fermion determinant. Assuming mass-degenerated quarks for simplicity 
reads 



SJT,0) = -2N t 



I ^^(^log[l + nWe-Vr] +trclog[l+a t We -^]Y (5 ) 



Here E p — \J p 2 + M 2 is the energy and M the constituent quark mass. As one can see 
the diagonal part of the Polyakov loop corresponds to consider chemical potentials for 
different color species. Large color gauge invariance is implemented by just averaging 
over group elements. 

The partition function character of the Polyakov loop can be appreciated considering 
a system with Nr dynamical quarks and an extra putative heavy quark (not antiquark) 
of mass inn at rest located at a fixed point and with fixed spin and colour a = 1, . . .N c . 
From Eq. © the change in the effective action is 

S q {Nf+l)-S q (N f ) = -2\og(l+Cl aa e- E »/ T )^-2e-' n »/ T Cl aa (6) 
yielding the partition function 



= l + (naa)2e- m "/ T = l + ^(tr c Cl)2e- m »/ T 



(7) 



after averaging over color degrees of freedom implied by DQ.. Thus we get 



1 

Wc 
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lim 



1 



m H - 



Z(Nf + K 
. Z(N f ) 



-1 



m H /T 



(8) 



If we use the HRG, Eq. (0Q), to evaluate the r.h.s. we reproduce the HRG result, Eq. © 
for the Polyakov loop, providing confidence on the assumed coupling to quarks. 

The action Sg{Q) would follow from gluodynamics but it is exponentially suppressed 
at low temperatures and the distribution of Q. locally coincides with the Haar measure. 
A convenient model to account for Polyakov loop correlations at different points and 
compatible with group integration at equal points 05ll is [@] 



(tr c a(jc)tr c a- 1 (y)) 5c = e- CT ^l/ r . 



(9) 



with o the string tension, including the correct screening of the color charge at large 
distances. This defines a correlation length and independent confinement domains with 
volume VfT = which describe the cross-over between deconfinement and chiral 

C° j 1 I 1 I 1 

symmetry restoration |j2l U22Ll34fl . 

Using these assumptions let us now consider the calculation of the partition function 
at low temperatures where to lowest order one gets contributions of just qq states 



which using Eq. © can be rewritten as 

io g z= (iN f f i d Y\! x d Y\r e - H ^™ ),T +■■■■ 

We recognize the classical partition function of a qq system with a Hamiltonian 



(ID 



H(x h p 1 ;x 2 ,p 2 ) = ^p 2 l +M 2 + ^/pl+M 2 + or l2 . 
Separating CM and relative motion, direct integration gives at low T values 



(12) 



logZ « (2N f ) 2 V J d\e-° r l J J S^e- E P' T 



(2N f ) 2 VV a 



MT 
2% 



-2M/T 



• (13) 



Similarly, for the Polyakov loop (heavy quark located at Xq) we get 

d 3 xd 3 p _ En/T ]_ 



2N 



. e -»Pi 



(tr c n(jc )tr c n l (x))s c + -- 



{lit? ' N c 
whence one can also rewrite the expression as 



(14) 



L(T) 



2N f [ d3xd3 P r -H(x,p)/T + . . . 



Nc 



{my 



(15) 



This corresponds to the classical partition function of the one-quark Hamiltonian 

H(x,p) = y/p 2 +M 2 + or. (16) 

where after integration 

, x N f M 2 T [M\ 2N f (MT \ 3/2 MlT 

As it was shown in our previous works these rules provide a satisfactory phenomenolog- 
ical description of the chiral-deconfinement cross-over observed in lattice calculations. 
The previous lowest order approximations do indeed reproduce the more sophisticated 
results up to T ~ 0.75r c [HI]. 



QUARK-HADRON DUALITY AT FINITE TEMPERATURE 

While the previous formulation reproduces lattice results in a satisfactory manner at 
T c ps 250MeV, within a hadronic phase it should be possible to express all observables 
in terms of purely hadronic properties, a feature absent in the model. To improve on 
this from Eqs. (fTTj) and (fT2j) we apply standard quantization rules to the relativistic two- 
body quantum-mechanical problem, which in the CM system (p l +p 2 = 0) corresponds 
to solve a Salpeter equation for scalar particles 

(2y / p 2 + M 2 + or^ y n = M n y n . (18) 

In a more elaborated treatment a two-body Dirac equation should be obtained. Note 
that at any stage our approach is Pauli-principle preserving at the quark level. Eq. (fT2l) . 
or its Dirac version, describes the interaction of qq pairs to form mesons, each meson 
with (2A^) 2 -fold degeneracy. The crucial point here is to keep track of the number 
and labeling of states contributing to the sum in Eq. (fTTT) . Thus after quantization and 
implementation of relativistic invariance, we get 

Here x,p represent the former CM coordinates of the qq pair. Note that the result holds 
even if quark masses are not degenerated. 

The expression obtained nicely reproduces the first bosonic term in the expansion 
of the RGM in Eq. (OQ). We are assuming that the pion is also contained in the sum 
although the dynamics leading to such a state is not literally determined by the specific 
Hamiltonian in Eq. (fT2~l) . The extension of the previous result to multiquark mesonic and 
baryonic states will be presented elsewhere Q36G . 

The argument can be extended to the Polyakov loop by re-quantization of the heavy- 
light Hamiltonian Eq. (TT6l) or its Dirac counterpart yieldin g th e eigenvalues A a and 
reproducing indeed the HRG form, Eq. (O, discussed in Ref. UM as required by quark- 
hadron duality. This constitutes our main insight which suggests re-analyzing these 
models after incorporating the connection to the HRG pal . 



The original Polyakov loop models suggested that Lj and Z (and hence (qq)r) are 
closely intertwined through their exponential suppression at low temperatures, con- 
trolled by the constituent quark mass (Eqs. (PT3l) and (fTTl)). However, in the hadronized 
version, the suppression depends on two not directly related hadron masses, namely, the 
pion for Z and the lightest meson with a heavy quark for the Polyakov loop. It is also 
noteworthy that the prefactors, namely, the powers of T present in the original model for 
both Z and L(T), are removed by the quantization. The origin of these prefactors was 
the relative motion in the qq pair, which produces a continuous spectrum in the classical 
case, but yields a discrete spectrum after quantization, typical of quantum bound states. 
Finally, the independence between deconfinement and chiral symmetry restoration cor- 
responding to crossed correlator (qqixe lA °l T )j ~ dL/dm q ~ becomes evident at low 
T due to the weak dependence of the Polyakov loop and thus of the static spectrum A a 
on the current quark mass. These features are also seen in lattice calculations 
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